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ABSTRACT 


Earthquake  source  theories  of  Haskell,  Brune,  and 
Savage  have  been  drawn  upon  to  develop  a description 
of  an  earthquake  as  a major  slip  accompanied  by  many 
smaller  tensional  and  slip  events.  We  find  natural 
explanations  of  several  previously  unexplained  obser- 
vations, such  as: 


Robustness  of  the  Ms:m^  discriminant, 

P corner  frequency  higher  than  S corner 
frequency 

High-frequency  P/S  amplitude  ratio  higher 
than  previous  theories  predict, 

Increase  of  complexity  as  a function  of 
third  moment, 


Small  m,  relative  to  M for  transform 
b s 


faults . 


The  theory  predicts  that  M^m^  populations  of  earth- 
quakes and  explosions  will  not  converge  at  small 
magnitudes . 
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INTRODUCTION 


Haskell  (1964)  developed  a theory  for  the  spectral 
distribution  of  the  tclcscismic  compressional  and  shear 
radiation  from  long  thin  strike-slip  or  tcnsional  earth- 
quakes in  a homogeneous  infinite  elastic  space.  He 
noted  that  the  strike-slip  earthquakes  had  a higher 
S/P  (shear-to-compressional)  energy  ratio  than  was 
commonly  observed,  and  speculated  that  most  earthquakes 
have  a small  component  of  tcnsional  faulting,  ile 
suggested  that  since  tensional  faults  have  a lower  S/P 
ratio  than  do  strike-slip  earthquakes,  a small  portion 
of  tcnsional  cracking  could  dramatically  alter  the  S/I5 
ratio.  Haskell  also  had  some  difficulty  in  matching 
the  existing  evidence  on  the  ratio  of  short-period  to 
long-period  energy,  and  speculated  that  the  fault 
surface  might  be  "rough".  This  possibility  he  modeled 
in  two  ways:  first,  by  allowing  the  ramp  displacement 

function  to  be  modulated  by  a sine  wave,  and  second, 
by  assuming  that  for  high  frequencies  the  fault  acts 
as  the  sum  of  several  small  faults.  Haskell's  (1969) 
paper  developed  techniques  for  calculating  the  near- 
field  displacement  from  slip  or  tensional  earthquakes. 

In  Haskell  (1966)  he  extended  the  idea  of  a 
modulated  ramp  displacement  function  to  a stochastic 
displacement  function.  The  ensemble  of  displacement 
functions  was  characterized  by  the  average  autocorre- 
lation of  the  ensemble.  This  average  autocorrelation 
led  then  to  an  average  spectrum.  Haskell  chose  the 
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autocorrelation  function  to  be  similar  to  that  of  the 
sine-wave  modulated  ramp  function  in  his  1964  paper, 
except  that  it  had  no  periodic  factor,  and  satisfied 
an  integral  constraint  resulting  from  the  fact  that 
the  earthquake  must  begin  and  end  in  a state  of  static 
equilibrium. 


The  displacement  spectrum  resulting  from  this 


autocorrelation  behaved  at  high  frequencies  like  w 


-3 


(the  spectrum  resulting  from  the  modulated  ramp  func- 


tion was  asymptotic  as  to  ).  Aki  (1967)  showed  that 

- 3 

an  co  asymptote  leads  to  predictions  for  the  shape 


of  the  Ms-m^  curve  in  disagreement  with  observation, 

and  chose  a different  ensemble  autocorrelation  func- 

_ o 

tion  which  yielded  an  asymptotic  slope  of  w „ (Aki's 
particular  choice  of  an  autocorrelation  does  not 
satisfy  Haskell's  criterion  that  the  earthquake  begin 
and  end  in  a state  of  static  equilibrium.  However, 
a different  autocorrelation  function,  which  would 

satisfy  the  equilibrium  criteria,  could  be  selected 

-2 


which  would  still  result  in  an  oj  asymptotic  slope.) 

Savage  (1966)  extended  Haskell's  1964  non- 
statistical  theory  to  handle  the  case  of  a dislocation 
nucleating  at  a point  in  the  fault  plane  and  spreading 
circularly  to  the  boundaries  of  a elliptical  fault. 
(Haskell  (1964)  had  assumed  that  the  fault  started  on 
a line  across  the  width  of  a narrow  fault  and  propa- 
gated down  the  length  of  the  fault).  Savage  showed 
that  while  the  simple  ramp  displacement  function 
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yielded  a displacement  spectrum  asymptotic  as  to  in 

Haskell's  initial  line  dislocation  model;  it  resulted 

in  an  to  asymptote  when  the  more  realistic  initial 

point  dislocation  was  assumed.  This,  he  showed,  was 

due  to  the  fact  that  the  radiating  surface  grew  quad- 

- 3 

ratically  with  time  instead  of  linearly.  This  to 
asymptote  was  corroborated  by  Molnar,  et.  al,  (1973) 
for  cases  in  which  the  fault  rupture  velocity  is  less 
than  the  shear  velocity. 

Brune  (1970)  presented  a shear  wave  source  theory 

which  brought  together  in  a physically  reasonable  way 

several  important  concepts.  He  made  plausible  the 

idea  that  until  influences  from  the  ends  of  the  fault 

were  felt,  the  fault  displacement  after  cracking  would 

be  a ramp  function  proportional  to  2-^-,  where  o is  the 

stress  and  y is  the  shear  modulus.  Assuming  then,  in 

effect,  that  the  entire  fault  began  to  radiate  at 

once,  the  characteristic  period  for  the  radiated  shear 

energy  became  L/S,  where  L is  the  fault  length  and  3 

is  the  shear  velocity.  Brune's  shear  wave  displace- 

_ 2 

ment  spectrum  decayed  asymptotically  as  w . This 
slope  was  not  strictly  speaking  derived,  but  simply 
followed  from  his  assumed  form  for  the  far-field  solu- 
tion which  he  chose  to  be  in  agreement  with  the  result 
for  an  instantaneous  stress  pulse  applied  to  the  in- 
terior of  a spherical  surface  (see  Bullen,  1963). 

This,  of  course,  is  not  necessarily  a good  model  of  a 
growing  displacement  nucleating  at  a point  on  the 
surface  of  a plane  or  spherical  surface.  (The  fact 


that  Bullen  [1963]  and  Randall  [1966]  found  corner  fre- 
quencies proportional  to  L / 3 for  shear  and  L/a  for  com- 
pressional  motion  reflects  the  fact  that  their  source 
was  assumed  to  act  instantaneously,  and  not  the  fact 
that  the  source  was  distributed  over  a volume  instead 
of  a surface,  as  suggested  by  Molnar  et  al.  [1973]). 

In  his  1970  paper  Brune  introduced  the  concept  of 
partial  stress  drop,  e,  and  showed  that  it  could  lead 
to  a spectral  slope  of  -1  somewhere  between  zero  fre- 
quency wher-1  the  displacement  spectrum  is  flat,  and 
asymptotic  frequencies  where  an  exponent  of  -2  obtains. 
Savage  (1972)  and  others  have  pointed  out  that  for  a 
long  thin  fault  one  must  expect  e<  W/L  where  W is  the 
fault  width  due  to  binding  of  the  fault  at  the  sides. 

Contemporaneous  with  his  theoretical  paper,  Brune 
and  a number  of  co-workers  analyzed  signals  from  many 
earthquakes  and  tried  to  interpret  them  in  its  light: 
see,  for  example,  Trifunac  and  Brune  (1970),  Wyss  (1970), 
Wyss  and  Brune  (1971),  Wyss  and  Hanks  (1972),  Hanks 
and  Wyss  (1972)  , Molnar  and  Wyss  (1972)  , Wyss  and  Molnar 
(1972).  One  of  the  most  interesting  findings  of  these 
analyses  is  that  the  P wave  corner  frequency  is  higher 
than  the  S wave  corner  frequency.  Hanks  and  Wyss  (1972) 
attempted  to  explain  this  fact  by  replacing  3 in  Brune’s 
(1970)  theory  by  a,  the  compressional  wave  velocity, 
and  appealing  to  the  concept  of  characteristic  wave- 
length. They  pointed  out,  however,  that  there  was  no 
theoretical  foundation  for  these  procedures.  Indeed 


Savage  (1972)  pointed  out  that  the  more  precise  theories 
of  Haskell  (1964)  and  Savage  (1966)  predict  approxi- 
mately equal  corner  frequencies  for  P and  S waves,  due 
to  the  fact  that  the  corner  frequencies  are  dominated 
by  the  duration  of  faulting,  so  long  as  that  time  is 
greater  than  L/B,  and  not  by  the  length  of  the  fault 
per  se. 

Molnar,  Tucker,  and  Brune  (1973)  have  reported 
that  the  P wave  corner  frequency  is  higher  than  that  of 
the  S wave  corner  for  144  aftershocks  of  the  San 
Fernando  February  9,  1971,  earthquake.  They  also  show 
by  an  exact  solution,  following  Savage  (1966),  that 
the  P-corner  is  higher  than  the  S-corner  if  the  strain 
is  released  instantaneously  over  the  entire  surface  of 
a circular  disk.  On  the  other  hand,  if  the  dislocation 
propagates  at  about  0.5  times  the  S wave  velocity  then 
they  find  that  the  two  corner  frequencies  are  about 
equal.  This  result  is  in  agreement  with  the  studies 
of  Haskell  (1964)  and  Savage  (1966,  1972)  . 

At  intermediate  dislocation  velocities  the  results 
of  Savage  (1974)  show  the  importance  of  the  operational 
definition  of  corner  frequency.  Savage  (1974)  shows 
that  if  the  corner  frequency  is  defined  as  the  inter- 
section of  the  high  frequency  displacement  asymptote 
with  the  horizontal  low- frequency  asymptote,  then  if 
v,  the  dislocation  velocity,  is  close  to  the  shear  wave 
velocity  the  P wave  corner  will  be  substantially  higher 
than  the  S wave  corner.  On  the  other  hand,  examination 
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of  Savage's  Figures  1 and  2 shows  that  for  the  same 
cases  the  P wave  corner  will  he  on  the  order  of  only 
201  higher  than  the  S wave  corner  if  the  corners  are 
taken  as  the  6dB  (1/2  amplitude)  points  on  the  dis- 
placement amplitude  spectrum.  The  identical  conclu- 
sions can  be  obtained  by  examination  of  Figure  6 in 
Molnar  et.  al.  (1973).  In  fact,  one  sees  there  the 
special  dangers  of  the  asymptotic  definition  of  corner 
frequencies  in  that  for  v=ot,  where  a is  the  compres- 
sional  wave  velocity,  the  intersection  definition  gives 
f <fg  for  0=85°,  an  observation  point  near  the  plane 
of  the  fault.  This  result  is  easily  verified  from 
Savage's  (1974)  figures,  and  one  may  see  also  that  for 
the  6dB  corners,  f >fg  . These  results  suggest  that 
the  6dB  corner  determination  is  more  reliable  than 
determination  by  the  intersection  of  asymptotes.  It 
is  probably  just  as  "objective"  as  the  intersection  of 
asymptotes  method,  and  is  certainly  not  as  severely 
affected  by  incorrect  values  of  Q,  a serious  problem 
with  the  latter  method  as  pointed  out  by  Thatcher  and 
Hanks  (1973). 


We  show  in  this  paper  that  another  explanation  of 


the  fp>fs  paradox  may  be  found  by  further  development 
of  Haskell's  (1964,  1966,  1972)  ideas  that  a substan- 


tial amount  of  tensional  faulting  accompanies  many 
strike-slip  or  dip-slip  earthquakes,  and  that  many 
earthquakes  are  accompanied  by  smaller  earthquakes. 
There  is  substantial  evidence  for  this  point  of  view 
in  the  literature.  Wyss  and  Brune  (1967)  concluded 
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that  the  28  March  1964  Alaskan  earthquake  ruptured  in 
a series  of  events;  initially  propagating  in  various 
azimuthal  directions  from  the  epicenter,  then  after  a 
period  of  about  44  seconds  propagating  steadily  to  the 
West.  They  also  concluded  that  other  large  earthquakes 
are  similar  in  character.  Trifunac  and  Brune  (1970) 
found  that  the  Imperial  Valley  earthquake  could  be 
regarded  as  a series  of  earthquakes  initiated  along  a 
strike-slip  fault  in  accordance  with  v/B~.5-.6.  Kamb 
et.  al.  (1971)  show  that  the  pattern  of  near-surface 
faulting  for  the  San  Fernando  earthquake  was  very 
complicated,  and  Jungels  and  Frazier  (1973)  conclude 
that  this  complexity  must  be  taken  into  account  in 
order  to  satisfactorily  interpret  the  static  displace- 
ment observations.  One  might  speculate  that  if  the 
near-surface  faulting  is  complicated,  so  oo  must  be 
the  deeper  faulting.  Jungels  and  Frazier  also  require 
changes  in  dip  for  the  fault  plane.  Bolt  (1972)  ana- 
lyzed the  Pacoima  Dam  strong  motion  records  for  the 
San  Fernando  earthquake  and  concluded  that  several 
"bursts"  of  high  frequency  energy  were  received  from 
remote  epicenters  on  the  fault  plane.  Mikumo  (1973) 
discusses  the  substantial  field  evidence  for  flexures 
in  the  San  Fernando  fault  plane  before  adopting  a 
plane  surface  for  his  theoretical  model.  Murray  (1973) 
has  concluded  that  the  main  Parkfield  earthquake 
consisted  of  two  separate  earthquakes  on  two  branches 
of  the  San  Andreas. 
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It  seems,  in  fact,  that  almost  every  carefully 
investigated  earthquake  reveals  first-order  departures 
from  a plant  slip  fault.  This,  of  course,  is  not  sur- 
prising giv^n  the  inhomogeneities  of  the  real  earth. 
Tensional  sub-earthquakes  must  also  be  expected. 
Although  a homogeneous  substance  may  fail  in  shear, 
one  must  expect  zones  of  weakness  in  an  inhomogeneous 
material  to  "pull  apart",  thus  introducing  a tensional 
component  to  an  earthquake.  Alternatively,  if  a fault 
plane  is  curved  and  strike-slip  movement  occurs,  one 
expects  some  separation  on  the  curves. 


The  above  remarks  are  not  meant  to  deny  that  the 
best  initial  approximation  to  almost  all  earthquakes 
is  slip  on  a plane.  This  is  strongly  implied  by  the 
numerous  successful  studies  of  long-period  radiation 
from  earthquakes  based  on  the  double -couple  model. 
However,  one  mi^ht  veasonably  expect  the  high-frequency 
displacement  spectrum  asymptotes  to  be  severely  af- 
fected by  small  tensional  sub-earthquakes. 

Another  puzzling  contradiction  between  theory  and 
observation  has  been  the  apparent  S/P  amplitude  ratio, 
as  observed  on  velocity  or  acceleration  seismograms , 
of  about  3/1,  whereas  theory,  e.g.  Haskell  (1964), 
predicts  a ratio  closer  to  10  or  20  to  1.  The  possi- 
bility that  observations  are  affected  by  unknown  vari- 
ations in  Q between  S and  P waves  somewhat  confuses 
the  issue.  The  most  reliable  data  is  presumably  that 
from  closest  in,  and  Trifunac  and  Brune  (1970)  show 
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strong-motion  seismograms  where  the  ratio  for  the  1940 
Imperial  Valley  subearthquakes  vary  between  10/1  and 
1/1  with  a mean  around  3/1.  The  first  six  minutes  of 
aftershocks  of  the  February  9,  1971  San  Fernando  Earth- 
quake as  recorded  at  the  Pacoima  Dam  strong  motion 
instrument  also  give  a mean  ratio  of  about  3/1,  as  we 
shall  see  in  later  sections  of  this  report.  As  pointed 
out  by  Haskell  (1964) , this  paradox  could  be  resolved 
if  much  of  the  energy  around  the  corner  frequency  came 
from  tensional  faulting  for  which  the  theoretical 
amplitude  ratio  is  about  3/1. 

Among  the  important  observations  which  a success- 
ful source  theory  should  address  are  those  which  have 
grown  out  of  the  efforts  to  discriminate  earthquakes 
from  underground  explosions.  The  amplitude  of  twenty- 
second  Rayleigh  waves  plotted  versus  the  amplitude  of 
one-second  compressional  waves  is  known  to  be  an  excel- 
lent discriminant;  see  for  example  Evernden  et.  al. 
(1969,  1971)  and  Marshall  and  Basham  (1972).  One  of 

Aki's  (1967)  chief  motivations  in  changing  from 

- 3 - 2 

Haskell's  (1966)  ui  model  to  an  w model  was  that 

the  latter  is  in  much  better  agreement  with  existing 

M -m,  data.  On  the  other  hand,  the  "spread"  in  observed 
s b 

^s'mb  pl°ts  5s  much  greater  than  can  be  accounted  for 
by  sampling,  measurement,  and  propagation  errors;  and 
therefore  any  single  line  in  M -m,  space,  such  as  the 
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oj  model  cannot  adequately  explain  the  data.  Examples 
of  especially  discordant  earthquakes  thac  look  like 
explosions  with  respect  to  this  discriminant  may  be 


* 


found  in  Landers  (1972)  and  Der  (197?).  Douglas  ev  al. 
(1973a)  have  shown  that  a substantial  amount  of  this 
discrepancy  may  be  explained  by  assuming  that  the 
earthquakes  are  dip-slip  with  a dip  of  45°.  This  fault- 
plane  orientation  has  a radiation  pattern  maximum  for 
the  teleseismic  P radiation.  Douglas  et  al.  and  Gilbert 
(1973)  have  also  shown  that  even  fcr  point  sources  the 
average  shallow  shearing  earthquake  will  separate  from 
shallow  explosions  by  approximately  0. 5-1.0  magnitude 
units.  A possible  physical  explanation  is  that  if  one 
first  matches  the  P wave  amplitudes,  then  for  the 
earthquake  there  is  ten  times  the  compressional  wave 
amplituda  emitted  as  shear  waves.  If  shear  waves  are 
about  as  efficiently  converted  to  Rayleigh  waves  at 
the  free  surface  as  are  the  compressional  waves  and  if 
the  earthquake  is  "shallow”  with  respect  to  the  Rayleigh 
wavelength,  then  one  expects  about  ten  times  the 
Rayleiph  wave  from  an  earthquake  as  from  an  explosion 
of  equal  m^ . The  same  argument  then  suggests  that  for 
tensional  earthquakes  one  would  have  only  about  three 
times  the  Rayleigh  wave  amplitude. 

Anglin  (1971)  plotted  complexity  versus  third 
moment  of  the  observed  spectrum  for  earthquakes  and 
explosions  and  found  good  separation.  In  general 
there  is  no  physical  explanation  for  the  separation, 
and  until  one  is  presented,  the  weaknesses  and  limita- 
tions of  the  discriminant  cannot  be  convincingly  dis- 
cussed. In  this  report  we  are  able  to  give  a plausible 
explanation  of  Anglin's  resuits  by  assuming  that  the 
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more  complex  the  earthquake,  the  more  small  earthquakes 
there  are  accompanying  it.  (A  complicating  point  in 
this  discussion  is  the  idea  developed  by  Douglas  et  al. 
(1973b)  that  complexity  may  be  due  to  arrivals  along 
separate  low  and  high-Q  paths.) 

In  this  report  we  find  that  the  high-frequency 
asymptotes  of  earthquakes  can  vary  between  w for 
simple  earthquakes  to  w ^ , where  the  exponent  is 

restrained  only  by  the  requirement  for  finite  radiated 
energy  (von  Seggern  and  Blandford  [1972]  present  evi- 
dence that  the  asymptote  for  explosions  is  id  ) .,  We 
find  that  in  the  case  of  complicated  earthquakes  mixing 
shear  and  tensional  sub -earthquakes , any  monotonically 
decreasing  displacement  spectral  shape  can  be  attained. 
Savage  (1972)  has  made  much  the  same  point,  saying 
"the  spectrum  of  the  incoherent  radiation  could  mask 
the  cjo  ^ asymptotic  behavior".  This  implies  that  short- 
period  discriminants  may  be  very  unreliable  in  practice. 
However,  with  knowledge  of  the  underlying  mechanisms 
it  may  be  possible  to  suitably  constrain  the  combina- 
tions of  discriminants,  regionalize  the  earthquakes, 
and  make  short  period  discriminants  predictable  and 
re  1 iab le . 

In  the  following  sections,  we  develop  and  apply 
the  theory,  make  some  comparisons  with  previous 
theoretical  results,  and  compare  some  results  with 
observation . 


i 
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THEORY 


Haskell  (1964)  gives  the  following  expressions 
for  the  teleseismic  radiation  from  a slip  fault  in  a 
homogeneous  full  space: 


ea(<u)  = (pw2L2Do/ 2tt3)  (B/a) 


G(u>)  |2Bi(u) 


e^Coj)  = (pw2L2D2/2ttB)  ) G(u))  | 2B2  (a)) 


(1) 


where  is  the  total  radiated  compressional  energy- 


density  at  frequency  ai,  a is  the  compressional  velocity, 
a)  is  the  circular  frequency,  p is  the  material  density, 
w is  the  fault  width,  L the  fault  length,  D0  the 
(constant)  displacement,  and  3 the  shear  velocity. 


E s satisfy  the  relation 
a, 3 1 


1 Z*00 

E Q = -* — f £ 
a, 3 2 it  Jo 


0 a, 3 


(u>)  da) 


where  E is  the  total  radiated  energy.  (Note:  To 

change  the  energy  spectrum  to  displacement  amplitude 

_ o 

spectrum,  multiply  by  io  and  take  the  square  root.) 

*G(o))  is  the  Fourier  transform  of  G(t)  where  G(t)  is 
the  displacement  time  history  at  a point  on  the  fault 
and  is  normalized  such  that  G=0  for  t<0  and  G->-l  as 
t-*-00.  Bi(o)),  illustrated  in  Figure  1 of  the  theoretical 

results  section,  is  asymptotic  to  0.27  at  low  frequencies 


- 2 

and  decays  like  u>  at  high  frequencies.  B2(u))  is 
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similar  but  has  a low  frequency  asymptote  of  0.40. 
ea  is  the  radiated  shear  energy. 

p 

These  are  Haskell's  relations  for  a long  thin 
fault  in  which  rupture  initiates  on  a line  across  the 
fault  and  for  which  the  direction  of  slip  is  down  the 
fault,  parallel  to  the  direction  of  rupture  propaga- 
tion. The  effects  of  the  velocity  of  rupture  propa- 
gation, v,  are  hidden  in  the  B functions  which  are 
defined  by,  for  example, 

BiO)  = |Trsin3ecos26X‘2sin2XC(de  (2) 

where 

Xq  = OjL (a  - cos 6) /2a  (3) 

and  a = a/v.  B2(w)  similarly  involves  b = B/v. 

Haskell’s  expressions  for  a fault  in  which  the 
displacement  is  in  the  fault  plane  perpendicular  to 
the  direction  of  rupture  propagation  a -e  identical 
except  that  B3  and  Bi*  replace  Bi  and  B2  respectively. 
We  shall  see  from  graphs  in  the  theoretical  results 
section  that  the  corresponding  functions  are  almost 
identical. 

For  tensional  faults  (in  which  the  displacement 
is  perpendicular  to  the  fault  plane),  the  expressions 
are  also  identical  except  that  B5  and  B6  replace  Bi 
and  B2;  and,  most  important,  the  expression  for  Ea 
has  been  multiplied  by  (a/3)4.  Since  the  ratio  B5/B6 
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is  almost  equal  to  the  ratio  B1/B2,  this  means  that 
the  displacement  S/P  ratio  for  tensional  earthquakes 
is  one-third  that  from  slip  earthquakes. 

Savage  (1966,  1974)  has  shown  that  the  corner 
frequencies  of  the  radiation  from  a fully  two-dimensional 
fault  are  determined  by  the  time  intervals  between  first 


motion  and  two  stopping  phases.  Each  corner  frequency 


contributes  a factor  00  to  the  far-field  displacement 
amplitude  spectrum.  In  turn,  we  may  see  by  an  argument 
much  like  Savage’s  (1966)  that  these  times  are  controlled 
by  the  two  principal  dimensions  of  the  fault.  Consider 
first  the  case  where  the  displacement  discontinuity 
(let  us  assume  a step  function)  nucleates  at  a point 
and  propagates  slowly  in  comparison  with  the  wave  velo- 
city. The  radiating  source  line  grows  linearly  in  time 
until  it  encounters  the  side  of  the  fault.  This  occa- 
sions the  first  stopping  phase,  since  the  radiating  line 
is  thereafter  constant  in  length  until  it  reaches  the  end 
of  the  fault  where  its  length  goes  to  zero  as  the 
wavefront  intersects  the  fault  edge.  This  generates 
the  final  stopping  phase.  Clearly  the  corner  frequencies 
will  be  proportional  to  v/W  and  v/L  for  both  P and  S 
waves . 


In  the  case  where  v>>a  and  B a similar  argument 
may  be  made.  Consider  those  receivers  not  near  the 
perpendicular  to  the  fault  plane  (as  Molnar  et  al. 

[ 19 7 3 J and  others  have  point  out,  these  receivers  fill 
the  dominant  portion  of  the  focal  sphere) . For  such 
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receivers  the  first  arrival  will  be  from  that  point 
on  the  edge  of  the  fault  nearest  the  receiver.  The 
line  on  the  fault  from  which  radiation  will  be  received 
will  grow  linearly  in  length  as  a function  of  time 
until  its  ends  encounter  the  sides  of  the  fault  at  a 
time  characteristically  about  W/ot  for  the  P wave  and 
W/$  for  the  S wave.  At  this  point  a stopping  phase 
is  created  as  the  generating  line  ceases  to  grow0  The 
line  thereafter  is  of  constant  length  until  it  encounters 
the  end  of  the  fault  at  times  about  L/a  and  L/3  for  P 
and  S respectively.  Thus,  as  assumed  by  Hanks  and 
Wyss  (1972)  and  others  and  as  shown  explicitly  by 
Molnar  et  al.  (1973)  and  Savage  (1974),  one  can  obtain 
different  P and  S corner  frequencies  if  v>“o  Whether 
this  is  a realistic  model  is  another  matter.  Note 
that  this  analysis  does  not  explicitly  call  on  dif- 
fraction theory  to  justify  the  location  of  the  corner 
frequencies.  With  further  geometrical  insight  one  can 
see  that  these  times  should  be  multiplied  by  sine, 
where  0 is  the  angle  between  the  perpendicular  to  the 
fault  and  the  direction  of  wave  propagation.  Other 
qualitative  aspects  of  the  radiation  can  be  deduced 
according  as  the  direction  of  wave  propagation  is 
almost  parallel  to  the  short  or  long  sides  of  the 
fault.  A rather  complete  picture  of  the  radiation  can 
be  built  up  from  physical  arguments  without  mathematics. 
See  Savage  (1966)  for  a detailed  discussion. 

We  can  now  see  that  Haskell*s  (1964)  analysis 
of  a long  thin  fault  assumes  in  effect  that  v/W  is  a 
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substantially  higher  frequency  than  those  under  obser- 

- 2 

vation.  We  can  also  see  that  contributing  to  the  to 
asymptotic  spectrum  obtained  by  Haskell  is  to  from 
the  B function  appropriate  to  L,  and  to  ^ from  his 
assumed  ramp  time  function.  It  is  clear,  as  pointed 
out  by  Savage  (1972)  and  confirmed  by  Molnar  et  al. 
(1973) , that  a ramp  time  function  in  connection  with  a 
realistic  fault  must  generate  an  asymptotic  spectrum 

_ 7 

proportional  to  10  . A minor  exception  occurs  in  the 

case  of  a circular  fault  where  the  displacement  nucle- 
ates at  the  center.  The  simultaneous  stopping  along 

the  whole  length  of  the  line  source  generates  a sharp 

..-2.5 

discontinuity  in  the  radiated  trace  resulting  in  to 
Savage  (1974,  1966).  Something  so  perfectly  symetrical 
would  be  unlikely  to  occur  in  the  earth. 

Archambeau’s  (1968)  relaxation  theory  also  results 
in  a spectrum  asymptotic  as  to  (Archambeau,  personal 
communication;  Minster  (1973).  However,  the  relation 
of  this  result  to  our  displacement  models  is  unclear. 

In  this  report  we  modify  Haskell’s  solutions  to 
take  account  of  a finite  value  for  W,  by  multiplying 
equations  (1)  by  an  additional  factor  and  B£  as 
appropriate;  except  that  L in  (3)  is  replaced  by  W. 

ejco)  = (pW2 L2Dq / 2rrB)  (3/ ot)  5|G(to)  | zB^(to)*B^ (to) 

n T w cn 

t (to)  = (pW2L2D§/2ttB)  |G(to)  ( 2Bz  (u>) -B2  (to) 

p 
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This  results  in  a solution  which  has  the  proper  quali- 
tative properties  such  as  generation  of  high-frequency 
snear  waves  when  v~6,  generating  waves  which  approxi- 
mate a shock,  corner  frequencies  as  appropriate  to  W 
and  L,  and  u)  ^ displacement  asymptotes.  In  view  of 
these  correct  qualitative  features  we  are  certain  that 
the  conclusions  reached  in  this  report  are  accurate. 

It  must  be  admitted,  however,  that  a more  satisfactory 
solution  would  be  to  use  solutions  such  as  those  pre- 
sented by  Molnar  et  al.  (1973)  and  Savage  (1974)  aver- 
aged over  the  focal  sphere  for  a suitable  non-circular 
fault  surface.  This  would  be,  however,  a difficult 
calculation.  We  have  therefore  fallen  back  to  an 
approximate  generalization  of  Haskell's  (1964)  solu- 
tion which  he,  of  course,  had  already  integrated  over 
the  focal  sphere.  We  might  remark  that  one  restriction 
in  Haskell's  analysis  results  from  the  fact  that  he 
represented  the  B integrals  as  the  sum  of  two  separate 
integrals,  each  of  which  diverges  if  v>B  but  whose  sum 
is  finite,  as  are  the  basic  integrals  (2).  Thus  in 
this  report  we  are  able  to  give  no  results  for  v _>B . 

Returning  to  equation  (1")  we  need  an  expression 
for  the  function  G.  In  this  we  follow  Brune  (1970) 
and  write : 

G(t)  = D ( t) / D (5) 


D = 0 , t<0 


0<t<£ " 
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t>e'  (6c) 


X = — 


L 

6 * 


e = 


eL 

21 


D = 


_ eagx 


eaL 

U 


" ^ (7) 


G(u) 


2 _ U) 


£ 2 T 2 (i)2  +T  “ 2 


(2-2e)  (l-cos(e'w))+e2)  . (8) 


Here  we  see  expressed  Brur.e’s  idea  that  a stress 
pulse  a is  applied  to  the  surface  of  the  fault  in  a 
medium  with  shear  modulus  p,  and  the  displacement  grows 
linearly  with  time  as  a ramp  function  until  approxi- 
mately the  time  t,  when  elastic  waves  return  from  the 
edges  of  the  fault  and  bring  the  motion  to  a gradual 
halt.  These  ideas  have  received  experimental  confir- 
mation in  work  by  Johnson,  Wu,  and  Scholz  (1973)  and 
by  Brune  (1973).  If  there  is  partial  stress  drop, 
then  at  a time  e"  a reverse  stress  is  applied  whose 
resulting  motion  also  comes  to  a completion  in  time  x 
and  reduces  tfte  final  displacement  from  (a/  )L  to 
e(a/p)L.  Savage  (1972)  and  others  have  pointed  out 
that  in  a long  narrow  fault  one  would  expect  ezW/L, 
since  the  absence  of  slip  at  the  edges  would  prevent 
a complete  dislocation.  In  such  a case  one  might 
better  choose  x=W/g  and  do  away  with  the  second  term 
of  equation  (6c)  unless  e<W/g,  The  expression  for  e' 
is  derived  by  assuming  that  the  reverse  stress  is  applied 
when  the  displacement  has  reached  1/2  its  final  value. 
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Our  expression  for  x is  in  agreement  with  that  of 
Brune  (1970),  but  the  expression  for  e'  is  different 
by  a factor  of  slightly  less  than  2 if  one  interprets 
the  radius  of  Brune* s circular  fault  as  equal  to  L/2. 
Such  a disagreement  is  easily  understandable  in  view 
of  our  slightly  different  formulation  of  the  problem, 
and  an  exact  answer  is  undoubtably  meaningless  since 
in  both  cases  the  motion  on  the  fault  is  prescribed 
and  is  not  a solution  of  the  dynamical  equations.  It 
seems  doubtful  that  it  is  possible  to  guess  the  dyna- 
mical motion  on  the  fault  to  better  than  a factor  of 
2.  We  shall  see  in  later  sections  that  an  error  in 
e'  of  a factor  of  2 would  have  no  significant  effect 
on  the  important  results  of  this  report. 

Equation  (8)  is  the  modulus  squared  of  the  Fourier 
transform  of  G(t).  The  expression  within  curly  braces 

is  the  square  of  Brune's  F(e).  We  see  that  the  asymp- 

••  2 o - ? 

totic  behavior  of  |G|  is  u) . Together  with  the  u> 

_ 4 

of  B this  yields  w for  e (with  two  B functions)  and 

- 3 a 

to  for  the  displacement. 

To  introduce  the  ideas  of  "subearthquakes"  devel- 
oped in  the  Introduction,  we  simply  sum  together  the 
energy  from  all  the  earthquakes  which  are  considered 
to  make  a contribution.  Here  we  are  assuming  that  the 
subearthquake  radiation  is  incoherent  with  respect  to 
all  other  earthquakes.  (For  identical  fault  orienta- 
tions and  slips  this  would  not  be  true  for  the  lowest 
frequencies.  We  assume  that  the  total  effect  of  such 
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subearthquakes  is  negligible  at  low  frequencies.)  The 

displacement  is  then  calculated  by  multiplying  the  sum 

- 2 

spectrum  by  u)  and  taking  the  square  root. 

Because  of  success  by  many  workers  in  analyzing 
the  long-period  radiation  from  earthquakes  in  terms  of 
the  simple  double-couple  model,  we  allow  the  moment  of 
the  earthquake  to  be  perturbed  by  subearthquakes  by 
substantially  less  than  a factor  of  2.  As  we  shall  see, 
this  places  practically  no  restrictions  on  the  shape 
and  character  of  the  higher  frequency  event  spectra. 

We  assume  that  the  subearthquakes  are  composed  in  eoual 
parts  of  slip  in  the  direction  of  their  own  fault  proga- 
gation,  slip  in  the  plane  perpendicular  to  the  direction 
of  fault  propagation,  and  tensional  motion  perpendicular 
to  the  fault  plane.  Each  of  the  three  subearthquake 
types  is  assumed  to  contribute  equally  to  the  submoment. 

In  principle  each  subearthquake  can  have  subearth- 
quakes of  its  own.  We  investigate  this  infinite 
regression  to  depths  up  to  n=6.  For  lack  of  a better 
assumption,  or  suitable  data,  we  assume  that  each  sub- 
level  entering  makes  an  equal  contribution  to  the  total 
submoment.  This  assumption  will  be  seen  to  require 
increasingly  many  earthquakes  at  each  sublevel. 

There  is  of  course  a distribution  of  the  size  of 
subearthquakes  ranging  from  the  size  of  the  main  earth- 
quake down  to  millimeters.  A sample  of  such  a distri- 
bution is,  however,  unknown  from  observation  or  theory 
and  so  we  assume  in  this  report  that  it  is  given  by  a 
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delta  function  at  some  fractional  length  (r  for  the 
n'th  level)  smaller  than  the  main  earthquake.  We 
assume  also  the  same  shape  and  the  same  values  for  a, 
y,  f, , a , v and  e.  Each  of  these  assumptions  could,  of 
course,  be  relaxed  if  data  or  a reasonable  hypothesis 
suggested  otherwise.  We  choose  them  only  because  they 
are  about  the  simplest  possible  assumptions.  Physically, 
they  amount  to  Aki’s  (1967)  similarity  hypothesis, 
although  we  make  the  hypothesis  here  for  subearthquakes 
rather  than  a set  of  separate  earthquakes.  Also,  each 
individual  earthquakes 1 s asymptotic  spectrum  decays  as 

_ 7 _ 2 

a)  instead  of  oo  as  in  Aki’s  theory. 

Given  the  assumptions  discussed  in  the  two  para- 
graphs above,  if  the  additional  submoment  is  M&,  and 
if  there  are  a total  of  N sublevels,  then  the  number 
of  earthquakes  of  a given  type  (two  slips,  one  ten- 
sional)  on  sublevel  n is  given  by: 

r ! 2 


Na  ’ 


eaL  w/y 


1 - 3n 

r 


We  have  explicitly  assumed  in  deriving  this  equa- 
tion that  the  low-frequency  "moment"  radiation  from  the 
subearthquakes  is  incoherent.  Differing  fault  plane 
orientations  could  easily  cause  this  even  if  the  faults 
were  all  slip  or  all  tensional. 

For  explosions  we  base  our  analysis  on  the  HARDHAT 
granite  data  of  Werth  and  Herbst  (1963)  as  analyzed  by 
von  Seggern  and  BLandford  (1972).  von  Seggern  and 
Blandford  followed  the  general  analysis  techniques  of 
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Haskell  (1967),  but  where  Haskell  only  allowed  the 
reduced  displacement  potential  to  imply  a step  in  the 
derivative  of  acceleration,  von  Seggern  and  Blandford 
allowed  a step  in  velocity.  This  assumption  seemed  to 
be  in  better  agreement  with  the  near-field  data  and 
resulted  in  an  asymptotic  displacement  spectrum  propor- 
tional to  a)  2 instead  of  of4  as  found  by  Haskell  (1967). 
This  scaling  also  allowed  von  Seggern  and  Blandford  to 
fit  teleseismic  data  from  Longshot,  Milrow,  and  Cannikin 
more  accurately. 

An  apparent  deficiency  of  Haskell's  (1967)  w"4 
explosion  model  which  prompted  von  Seggern  and  Bland- 
ford' s work  was  the  fact  that  at  high  frequencies 
explosions  gave  a decreasing  amplitude  spectrum  as  a 
function  of  yield.  Thus  for  large  enough  yield,  cu^e 
root  scaling  would  give  decreasing  m,  values.  A simi- 

*7  D 

lar  phenomenon  occurs  in  the  w earthquake  spectra  in 
this  report.  One  obtains  a constant  m^  as  fault  length 
increases.  At  a lew  magnitude  the  N^-m^  curve  turns 
up  sharply,  in  disagreement  with  observation.  Haskell's 
(1966)  a)  statistical  earthquake  model  implied  the 

same  disagreement  and  this  prompted  Aki  (1967)  to  pro- 
_ 2 r 
pose  his  a)  model.  In  this  report  we  satisfy  observa- 
tion by  allowing  subearthquakes. 

The  expression  for  the  explosion  far-field  energy 
spectrum  from  von  Seggern  and  Blandford  (1972)  is: 
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e*(u})  = ( 8TT p ‘P 2 (°°) u)2/a)  IC^B^/k) — +lj_ 
“ [ (w/k) 2+l] 3 


n°o  = t0 

o 


k = k 


Y 

7" 


-1/3 


(10) 


Here  ^(o®)  is  the  long-time  asymptotic  reduced 

displacement  potential,  B is  a dimensionless  constant, 

and  k is  a characteristic  frequency.  Cube-root  scaling 

implies  the  other  relations  in  (10),  where  Y and  k 

o o 

are  observed  values  from  a scaling  event.  Using  Werth 
and  Herbst's  (1963)  data  for  HARDHAT,  von  Seggern  and 
Blandford  find  B = 2.04,  V = 2.5  x 109cm3,  Y = 5 


kilotons,  k = 16.8  sec-1. 


From  Haskell  (1967)  we  have  for  the  explosion 
displacement  spectrum  at  radius  R, 


UX  = (2pa47rR2u)2) _1/2  [^(w)]1/2. 


(ID 


The  corresponding  expressions  for  earthquakes 


are : 


= (2pa47TR2u)2)’1/2  [ea(u>)]1/2 


(12) 


U®  = (2  p(3  4ttR2u)2)  " 1/2  [ee(U)]1/2( 
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THEORETICAL  RESULTS 

In  Figure  1 we  see  B1  and  B2  presented  as  a func- 
tion of  p = toL/a  or  toL/3  as  appropriate  „ We  choose 
v = 0.7753  and  a = /33  so  that  these  are  identical  to 
the  functions  graphically  displayed  by  Haskell  (1964)  . 

In  Figure  2 we  compare  the  B functions  for  the 
three  types  of  fault  discussed  by  Haskell  (1964)  . We 
see  that  they  are  quite  similar  and  could  probably  be 
replaced  by  one  function  for  P and  another  function 
for  S for  most  qualitative  and  quantitative  purposes. 

It  might  be  suggested  that  in  equations  (1“*)  we 

2 2 

should  use  B^B^  and  B2B4  instead  of  B^  and  B2,  since 

for  a fault  nucleating  at  a point,  half  the  fault 

slips  in  the  direction  of  propagation,  and  half  at 

« 

right  angles.  However,  inspection  of  Figure  2,  in 
connection  with  other  results  to  be  discussed,  will 
show  that  the  net  effect  would  be  negligible. 

Figure  3 illustrates  the  effect  of  increasing  or 

decreasing  the  velocity  of  rupture.  As  mentioned  in 

the  Theory  section,  Haskell's  analytical  approach  makes 

it  impossible  to  calculate  solutions  for  v _>  3.  One 

easily  sees  that  lower  velocities  reduce  the  radiated 

energy  and  decrease  the  S/P  ratio.  Another  interesting 

-1 

feature  is  the  introduction  of  an  apparent  to  inter- 
mediate trend  for  v z 3«  This  probably  means  that  the 
low-frequency  shear  waves  have  formed  a shock  wave, 
but  the  higher  frequency  waves  are  still  not  in  phase. 
Similar  behavior  in  the  near-field  was  noted  by  Aki 
(1968) . Thus  we  have  another  cause  for  an  intermediate 
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trend  in  the  displacement  spectrum,  in  addition  to  fault 
width  and  partial  stress  drop.  Among  authors  who  have 
discussed  the  possible  velocities  of  fault  rupture  are 
Mansinha  (1964),  Aki  (1968),  Weertman  (1969),  Wyss  and 
Brune  (1967),  Murray  (1973),  and  Boore,  Aki,  and  Todd 
(1971) . 

Figures  4 and  5 illustrate  the  same  conclusions 
for  the  other  two  types  of  fault  except  that  the  inter- 
mediate frequency  trend  is  not  so  apparent . 

Next  we  carry  the  analysis  to  the  displacement 
spectrum  given  by  (11).  Throughout  this  report  the 
earthquake  parameters  given  in  Table  1 are  held  fixed. 
The  parameters  L,  v,  e,  n (the  number  of  levels  of  sub- 
earthquakes), r (the  ratio  of  lengths  of  successive 
subearthquakes),,  and  earthquake  type  vary,  and  are 
indicated  as  appropriate  on  each  Figure.  We  choose 
the  letters  SS,  D5  and  T to  indicate  the  main  earth- 
quake type;  respectively,  slip  in  direction  of  propa- 
gation, perpendicular  to  direction  of  propagation  but 
in  fault  plane,  and  perpendicular  to  the  fault  plane. 
Solutions  for  other  values  of  a may  be  obtained  by 
scaling  the  displacements  up  or  down  in  proportion  to 
(a/100  bars)  . 

Table  1 

Earthquake  Parameters  Fixed  in  this  Report 

8-2  - 3 

a=  100  bars  = 10  dyne  cm  a=8.0km/sec  p = 3.0gr  cm 

ii  ? 3 

11=3.0x10  dyne  cm  3 = 4.62km/sec  R=8.0xl0  km 

W=L 


In  Figure  6 we  have  plotted  the  displacement  spectra 
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for  P and  S,  and  their  ratio.  Note  that  the  corner  frequen- 
cies for  S and  P (indicated  by  filled  circles  on  this  and 
subsequent  Figures)  are  almost  equal.  Note  also  that  the 
P/S  ratio  is  about  0.13  decreasing  to  about  0.05  for  higher 
frequencies.  The  sharp  dip  around  1.7  Hz  is  probably 
not  realistic,  for  two  reasons:  first,  as  we  noted 

above,  use  of  B.^  instead  of  B.^  is  probably  more 
realistic  and  would  smooth  out  this  null;  second,  in 
the  real  earth  any  null  this  sharp  is  usually  smoothed 
out  by  effects  not  considered  here. 

Next  we  investigate  the  effect  of  subearthquakes 
by  adding  to  the  simple  type  SS  fault  just  discussed 
equal  numbers  of  fault  types  SS,  DS,  and  T with 
lengths  = L/3  = 1.05  km.,  in  sufficient  quantity  to 
increase  the  long  period  P displacement  radiation  by 
about  25%;  and  the  long-period  S radiation  by  about 
10%.  We  see  in  Figure  7 that  as  a result  w is  higher 
than  u>s  by  a factor  of  1.7,  and  the  high-frequency 
P/S  ratio  varies  between  0.17  for  low  frequencies  to 
0.28  for  high  frequencies.  Note  also  that  it  would 
be  extremely  difficult  to  determine  ws  by  use  of  the 
high-frequency  asymptote  since  the  spectrum  has  a 
rather  uniform  curvature  at  high  frequencies.  We  see 
also  in  Figure  7 the  definite  suggestion  of  an  inter- 
mediate slope.  This  last  effect  of  the  subearthquakes 
would  be  created  no  matter  what  their  type.  But  the 
shift  in  relative  corner  frequency  and  the  increase 
in  P/S  ratio  at  high  frequencies  are  due  solely  to  the 
presence  of  1/3  tensional  subearti  quakes . Since  their 
P energy  is  comparatively  high,  and  since  earthquakes 
of  any  size  of  a given  type  have  equal  high-frequency 
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displacements  (because  of  the  lo  asymptotes)  , almost 
all  the  high-frequency  P wave  energy  comes  from  the 
several  small  tensional  subearthquakes.  Therefore  oj 
is  larger  than  whose  value  is  not  so  severely 
affected  by  the  subearthquake  energy. 

Figure  8 shows  explicitly  how  simple  earthquakes 
of  different  sizes  have  identical  high-frequency  dis- 
placement spectra.  The  "periodicity"  in  the  P/S  ratio 
illustrates  the  fact  that  exactly  the  same  B functions 
are  shifted  left  and  right  along  the  frequency  scale 
to  produce  plots  for  different  size  earthquakes. 

Figure  9 shows  e/plicitly  the  P/S  ratio  for 
different  types  of  earthquakes.  We  see  that  at  all 
frequencies  the  tensional  earthquakes  he  2 about  3 
times  the  compress ional  amplitude  of  slip  earthquakes. 

Figure  10  adds  to  the  energy  spectra  of  Figure  9 
the  energy  spectra  from  equal  numbers  of  fault  types 
SS,  DS,  and  T with  lengths  L/3  = 1.05  km  in  sufficient 
quantity  to  increase  the  long-period  P displacement  by 
about  25%,  and  the  long-period  S radiation  by  about  10%, 
We  see  that  the  spectra  including  the  SS  and  DS  prin- 
cipal earthquakes  become  identical  to  each  other, 
the  same  as  we  have  seen  already  in  Figure  7.  The 
spectrum  including  the  principal  T earthquake  is  equal 
to  the  other  two  above  1.0  Hz,  but  below  1.0  Hz  has 
about  3 times  the  P wave  amplitude. 


In  Figure  11  we  investigate  the  effects  of  diffe- 
rent values  for  n.  Recall  from  the  Theory  section 
that  we  force  each  sub- level  to  make  an  equal  contri- 
bution to  the  total  moment,  but  do  not  allow  the  total 
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moment  to  change  as  n increases;  see  equation  (9).  As 
might  be  expected  we  are  able  to  move  the  P wave  corner 
frequency  arbitrarily  high.  It  seems  clear  also  that 
earthquakes  of  greatly  differing  size  from  a single 
region  might  have  the  same  corner  frequency  if  the 
smallest  subearthquakes  were  always  the  same  size, 
perhaps  being  controlled  by  the  lack  of  rock  joints 
below  a certain  scale. 


In  Figure  12  we  vary  the  total  moment  of  the 
subearthquakes  and  see  that  once  the  total  moment  rises 
above,  say  10%  of  the  original,  the  effect  on  corner 
frequency  is  established,  and  further  increases  do 
not  seem  to  change  the  situation  substantially.  It  is 
noticable,  however,  that  subearthquakes  that  create 
undetectable  amplitude  changes  at  frequencies  as  high 
as  0.4  Hz  can  have  an  important  influence  on  P/S  ratios 
around  2.0  Hz. 


In  Figure  13  we  vary  the  value  of  r while  main- 
taining the  total  moment  constant  and  using  only  n=l. 

As  might  be  expected,  we  see  that  the  smaller  earth- 
quakes have  the  more  profound  e “feet  on  corner  frequencie 
and  P/S  ratios  (of  course,  were  their  number  not  in- 
creased according  to  (9)  in  order  to  maintain  constant 
moment,  the  smaller  subearthquakes  would  have  a smaller 
effect)  . 


Figure  14  shows  the  effect  on  simple  earthquakes 
of  variation  in  v,  the  rupture  propagation  velocity. 

As  we  saw  in  connection  with  Fi'*  re  4,  an  intermediate 
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cj"1  spectral  t/end  with  frequency  emerges  at  the 
higher  velocities.  Note  also  that  the  large  velocities 
create  very  low  P/S  ratios  at  high  frequencies  - on  the 
order  of  0.01.  This  occurs  because  of  the  great  en- 
hancement of  the  S wave  spectra  in  shock  formation. 

In  Figure  16  we  introduce  variation  in  e for  a 
simple  earthquake  and  see  that  shift  in  moment  and  the 
appearance  of  an  intermediate  spectral  trend  with  fre- 
quency as  shown  by  Brune  (1970).  As  can  easily  be 
seen  from  the  structure  of  equation  (1'’),  changes  in 
e do  not  affect  the  P/S  ratio,  and  this  can  be  seen 
explicitly  in  Figure  16.  From  this  fact  one  would 
expect  that  while  the  P and  S corner  frequencies  might 
shift  with  e,  they  would  not  shift  relative  to  ono 
another.  We  see  from  Figure  16  that  this  is  in  fact 
true,  and  also  that  the  corner  frequencies  themselves 
change  by  less  than  201. 

Figure  17  illustrates  the  results  when  one  level 
of  subearthquakes  with  the  appropriate  values  of  e 
are  added  to  the  simple  earthquakes  of  Figure  16. 

We  might  repeat  at  this  point  that  solutions  for 
other  values  of  o or  p can  be  obtained  by  shifting 
lines  in  Figures  6-17  up  or  down  the  page,  since  these 
parameters  appear  only  as  multiplicative  constants; 
see  equations  (1",  5,  7,  8).  This  also  shows  explicitly 
that  earthquakes  of  greatly  different  magnitudes  may 
have  the  same  corner  frequency. 
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Next  we  discuss  the  theoretical  results  in  terms 
of  their  impact  on  the  Ms~m^  discriminant.  Rather 
than  try  at  this  point  in  the  report  to  relate  the 
displacement  values  calculated  at  .05  and  1.0  Hz  to  M 
and  m^  values,  we  shall  first  make  whatever  deductions 
are  possible  by  simply  assuming  proportionality  between 
log  displacement  and  magnitude.  We  shall,  however, 
slightly  anticipate  later  arguments  by  plotting  the 
root-mean-square  compres s ional  plus  shear  displacement 
for  the  0.05  Hz  displacement.  All  conclusions  reached 
in  this  report  would,  however,  have  been  identical  had 
we  used  only  the  shear  displacement,  since  the  root- 
mean-square  is  only  the  order  of  51  greater.  For  the 
1.0  Hz  displacement  we  use,  of  course,  only  the  com- 
pressional  wave  displacement. 


In  Figure  18,  we  see  the  M^-m^  relation  for  simple 
earthquakes.  For  the  smaller  fault  lengths  where  the 
corner  frequency  is  above  1.0  Hz,  the  spectrum  is  flat 


between  0.05  and  1.0  Hz  and  therefore  is  propor- 


tional to  m,  . Between  L = 1.0  and  3.16  km  the  corner 
b 


frequency  passes  through  1.0  Hz  and  the  displacement 
at  1.0  Hz  becomes  nearly  constant.  With  increasing  L 
the  0.05  Hz  displacement  continues  to  increase.  This, 


of  course,  results  in  a nearly  vertical  N^-m^  line. 


Around  L = 31.6  kilometers  the  corner  frequency  becomes 


lower  than  0.05  Hz,  and  Mg  no  longer  increases 


s ignif icantly . 


The  resulting  N^-m^  line  is  strikingly  different 
from  those  observed,  e.g.  Evernden  et  al.  (1969,  1971). 


I 
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As  mentioned  before,  such  a line  is  a necessary  result 
_ 3 

of  any  to  model.  This  disagreement  with  observation 


drove  Aki  (1967)  to  consider  to  models.  Let  us  note 
at  this  point,  however,  that  we  have  assumed  a = 100 
bars  for  all  calculations  in  this  report.  Thatcher 
(19  72),  Thatcher  and  Hanks  (19  73),  and  others  have 
documented  a values  ranging  roughly  from  0.2  to  200 
bars.  Hanks  (1974)  has  deduced  by  analysis  of  the 
Pacoima  Dam  seismograms  stress  on  the  order  of  350-1400 
bars  near  the  epicenter  of  the  San  Fernando  earthquake. 
Accordingly,  any  point  on  Figure  18  could  be  moved 
diagonally  as  much  as  3-4  magnitude  units.  It  can  be 
seen  that  this  corresponds  to  the  fact  that  variation 


in  a produces  a straight  line  in  the  ft  -f  (moment- 

oo 


corner- frequency)  diagrams  of  Hanks  and  Thatcher  (1972) 
and  Thatcher  and  Hanks  (1973),  Thus,  even  a region 
with  only  simple  earthquakes  could  have  a wide  range 
of  possible  Ms"m^  plots  if  there  were  a wide  range  in 
possible  values  for  a.  Note,  however,  that  none  of 
the  earthquakes  would  cross  the  limiting  line  which 
is  drawn  so  that  if  the  shear  and  compressional  spectra 
are  flat,  the  earthquake  will  lie  upon  it. 


In  Figures  19-22  we  show  the  effects  on  the  M -m, 

s b 


plot  of  2,  4,  or  6 levels  of  subearthquakes.  We  con- 
sider subearthquakes  for  which  r = 0.33  or  0.57,  and 
in  which  the  additional  moment  is  constrained  to  be 
equal  to  those  indicated  by  lines  3 or  5 in  Figure  12. 
These  two  levels  represent  percent  changes  in  the  P 
wave  long-period  radiation  of  25%  and  1%  respectively. 
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ive  see  that  above  the  limiting  line  almost  any  M -m 
curve  desired  may  be  achieved,  especially  when  one 
considers  possible  variations  in  a.  This  can  occur 
even  though  the  amplitude  of  the  P wave  long-period 
displacement  changes  by  less  than  1%.  Still,  we  see 
that  there  are  no  major  excursions  over  the  limiting 
Ms"mb  line-  The  few  small  violations  that  do  occur 
may  be  traced  to  the  influence  of  tensional  subearth- 
quakes. 

In  Figures  23-25  we  see  that  with  n=4  a variation 
of  £ and/or  v can  give  a rather  reasonable  M -m.  curve, 
and  with  variations  in  a might  be  adjusted  even^closer 
to  observation.  Obviously  we  have  embarrassingly  many 
variable  parameters,  at  our  disposal.  However,  in  none 
of  these  cases  is  there  a significant  excursion  over 
the  limiting  line. 

In  Figure  26  we  arrive  at  just  such  a significant 
excursion  by  assuming  that  the  basic  simple  earthquake 
is  tensional  in  character.  It  is  in  fact  easy  to  see 
that  the  limiting  line  for  tensional  earthquakes  is 
about  0.5  magnitude  units  to  the  right  of  (or  0.35  units 
perpendicular  to)  the  limiting  line  for  slip  earth- 
quakes. There  have,  however,  been  no  tensional  earth- 
quake fault-plane  solutions  reported  in  the  literature, 

so  that  the  existence  of  such  simple  earthquakes  is  in’ 
doubt. 

Finally,  in  Figures  27  and  28  we  see  the  explosion 
spectra  and  curves  resulting  from  equations  (10,11). 
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As  an  insight  into  the  next  section  about  comparison 
with  observation  we  have  placed  on  Figure  28  the 
limiting  M^-m^  line  for  slip  earthquakes.  A dashed 
reference  line  has  also  been  drawn  through  the  family 


of  explosion  points  to  show  that  Mg  is  not  strictly 


proportional  to  m^  even  at  intermediate  frequencies. 

In  fact,  we  notice  that  the  slope  is  less  than  1.0  for 
intermediate  yields,  a result  of  the  existence  of  the 
maximum  in  the  spectra  of  Figure  27.  Brune's  displace 
ment  function  (6)  results  in  a no  such  intermediate 
frequency  maxima  for  earthquakes.  If  one  were  to 
specify  a displacement  time  history  such  as  Haskell's 
(1964)  modulated  ramp  because  of  some  such  physical 
idea  as  overshooting  of  the  equilibrium  displacement, 
then  one  would  have  another  mechanism  for  obtaining 
earthquake  Ms_m^  points  on  the  lower  side  of  the 
limiting  line.  Archambeau's  (1968)  theory  apparently 
predicts  such  an  overshoot;  however,  it  has  not  been 
found  in  the  experimental  studies  of  Johnson  et  al. 
(1973)  and  Brune  (1973) . 
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COMPARISON  WITH  OBSERVATION 


Douglas  et  al.  (1973a)  have  calculated  theoretical 
values  for  explosion  and  also  assuming,  as  we 
have,  Werth  and  Herbst’s  (1963)  granite  reduced  dis- 
placement potential.  They  took  into  account  the 
distance  to  the  free  surface,  receiver  and  source 
region  layering,  absorption,  dispersion,  divergence  of 
the  P waves,  instrument  response,  and  conventional 
measurement  techniques  for  and  m^ . For  100  kilotons 
in  granite  at  a depth  such  that  pP  reinforces  P they 
obtained  mb  = 5.2,  Mg  =4.0.  If  we  subtract  0 . 2 mb 
units  as  an  approximation  to  the  mb  which  would  be 
obtained  for  a deep  shot  where  pP  did  not  reinforce  P, 
then  we  may  appropriately  calibrate  Figure  28  by 
setting  the  100  kt  shot  at  the  mb  = 5.0,  Mg  = 4.0 
point.  The  result  may  be  seen  in  Figure  29.  This  is, 
of  course,  a purely  theoretical  calibration,  and 
effects  such  as  spall,  anomalous  absorption,  etc., 
may  invalidate  it. 


To  gain  further  insight  into  the  nature  of  our 
solutions  we  superimpose  Figure  29  on  a figure  from 
Douglas  et  al.,  resulting  in  Figure  30.  The  para- 
meters for  the  earthquakes  plotted  in  Figure  30  are 
given  in  Table  2 taken  from  Douglas  et  al.  Some  of 
these  points  were  misplotted  in  their  original  figure, 
and  the  appropriate  corrections  have  been  made  in 
Figure  30.  We  see  that  some  explosion  points  fall  to 
the  left  of  our  granite  line.  This  may  be  traced  to 
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cancellation  of  P by  pP  from  explosions  buried  at  the 
minimum  containment  depth.  Thus  in  this  way  explosions 
can  be  made  "anomalous",  at  the  possible  risk  of  sub- 
sidence craters,  surface  leakage  of  radioactivity,  etc. 
Note  in  Figure  30  that  explosions  in  other  materials 
are  less  earthquake- like , although  the  corresponding 
reduced  displacement  potentials  may  not  be  as  reliable. 

The  earthquake  models  assumed  by  Douglas  et  al. 
were  typically  2 km  in  diameter;  the  assumed  displace- 
ment time  history  was  a step  nucleating  at  a point  and 
spreading  to  the  perimeter  of  an  ellipse  (basically 
Savage's  1966  model).  The  step  time  history  will 
cause  high  frequency  displacement  to  fall  off  propor- 
tional to  to'1,5  for  the  circular  fault  as  discussed 

7 

by  Savage  (1974)  instead  of  u as  in  our  model.  Thus 

the  2 km  diameter  faults  of  Douglas  et  al.  correspond 

most  closely  in  Figure  18  to  the  1 km  length  fault. 

Notice  that  this  is  just  at  the  break  point  where  the 

M -m,  line  for  simple  earthquakes  begins  to  depart 
s b 

dramatically  from  observation.  We  may  therefore  regard 
the  earthquake  calculations  of  Douglas  et  al.  to  be 
for  point  sources  with  respect  to 

Thus  the  fact  that  our  limiting  line  lies  below 
but  close  to  most  of  the  earthquakes  justifies  the 
implicit  assumption  we  discussed  in  the  Introduction. 
That  is,  for  point  sources  when  P wave  amplitudes  are 
matched,  earthquakes  will  have  a larger  Mg  because 
they  emit  an  equal  amount  of  long-period  P as  do 
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explosions,  plus  about  ten  times  the  P wave  amplitude 
as  shear  waves.  If  we  may  assume  the  absence  of  radi- 
ation pattern  effects,  and  that  S waves  convert  to 
Rayleigh  waves  as  well  as  do  P waves,  then  we  may 
anticipate  a full  magnitude  separation  in  the  M^-m^ 
lines. 

The  earthquakes  which  do  cross  the  limiting  line 
in  Figure  30  illustrate  the  weaknesses  in  the  assump- 
tion of  no  effects  due  to  radiation  pattern.  The 
triangles  correspond  to  earthquakes  at  a depth  such 
that  there  is  a null  in  the  20-second  Rayleigh  radia- 
tion. This  misclassification  would  be  overcome  by 
measurement  of  the  maximum  Rayleigh  motion  instead  of 
restricting  the  measurement  narrowly  to  the  20-second 
period.  The  low  (Mg-m^)  earthquake  at  m^  = 4.8, 

Mg  = 3.9  corresponds  to  a 43°  dip-slip  event  which  has 
a maximum  in  the  P radiation  for  teleseismic  distances. 
Such  events  would  also  be  especially  difficult  to 
discriminate  because  they  can  display  compres s ional 
first  motion  at  all  teleseismic  stations.  Since  rare- 
factions have  been  observed  at  distances  of  5°  from 
Der's  (1973)  low  (Mg-rn^)  events,  the  events  cannot  be 
purely  tensional  in  character.  This  leaves  the  possi- 
bility of  a complex  45°  dip-slip  event,  or  possibly 
some  unknown  path  or  depth  effect. 


In  Figure  31  we  have  superimposed  our  theoretical 
lines  on  the  data  of  hvernden  et  al.  (1971).  We  see 
that  most  of  the  explosions  lie  to  the  left  of  the 
theoretical  granite  line.  This  is  probably  because 


they  were  fired  at  the  minimum  containment,  depth, 
resulting  in  interference  between  pP  and  P.  The  fact 
that  the  largest  explosions  do  not  follow  the  granite 
line’s  curvature  is  probably  due  to  the  fact  that  m^ 
is  being  measured  at  frequencies  less  than  1.0  Hz, 
and  also  because  the  greater  confining  depth  has 
reduced  the  effect  of  the  pP-P  cancellation. 

The  limiting  line  for  earthquakes  seems  to  be  in 
good  agreement  vsith  the  observations.  This  suggests 
that  few  of  the  earthquakes  in  the  sample  were  domi- 
nated by  tensional  or  45°  dip-slip  faulting,  or  had 
displacement  time  histories  characterized  by  an  over- 
shoot of  the  equilibrium  position. 

Inspection  of  Figure  18  shows  that  there  is  a 
limit  of  about  4.4  Mg  units  between  the  limiting  line 
and  the  highest  possible  earthquake.  In  Figure  31, 
this  difference  is  only  about  1.0  magnitude  unit, 
suggesting  that  there  are  not  many  simple  earthquakes 
in  the  Southwest  United  States.  ^s_mb  diagrams  by 
other  authors,  e.g.  Marshall  and  Basham  (1972)  show  a 
spread  of  two  magnitude  units  for  Eurasian  events. 

Since  even  earthquakes  of  L=10  (see  Figure  18)  imply 
AMs^3,  there  are  apparently  few  simple  earthquakes  in 
Eurasia. 

The  success  of  this  comparison  with  observation 
suggests  that  the  proper  way  to  align  plot  of 

events  from  regions  where  there  have  been  no  explosions 
is,  after  station  corrections  have  been  determined  and 
applied,  to  draw  the  limiting  line  and  align  it  with 
some  standard  such  as  the  one  in  Figure  31. 
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Next  we  consider  the  questions  on  observations  of 
S/P  amplitude  ratios  and  their  connections  with  corner 
frequency.  Probably  the  best  data  available  consist 
of  the  Pacoima  Dam  records  of  the  first  six  minutes 
of  aftershocks  of  the  San  Fernando  earthquake.  Both 
P and  S waves  are  accurately  recorded  because  an  0.5° 
tilt  caused  by  the  major  earthquake  closed  the  instru- 
ment's trigger  and  caused  the  instrument  to  run  con- 
tinuously. In  addition,  the  epicenters  were  only  about 
10  kilometers  distant  and  10  kilometers  below  the 
surface,  so  that  effects  of  absorption  and  interaction 
with  the  free  surface  were  at  a minimum. 

The  acceleration  seismograms  are  shown  in  Figure  32, 
taken  from  Trifunac  (1972).  The  brackets  indicate  the 
P and  S wave  windows  used  by  Trifunac  in  his  accelero- 
gram spectral  analyses.  These  spectra  are  displayed 
in  Figure  33.  Even  a cursory  inspection  of  Figure  32 
shows  that  the  P/S  ratio  is  much  closer  to  0.3  than 
to  0.03,  as  predicted  by  the  theory  of  the  simple 
earthquake.  More  quantitatively,  measurement  of  the 
levels  of  the  high-frequency  asymptotes  in  Figure  33 
shows  that  the  median  ratio  is  0.37. 

Inspection  of  the  records  also  shows  that  the  P 
waves  are  of  higher  frequency  than  the  S waves,  and 
Figure  34  (from  Molnar  et  al.,  1973,  after  the  data 
in  Trifunac,  1972)  shows  this  explicitly.  In  the  con- 
text of  the  present  theory  these  results  imply  that 
these  San  Fernando  aftershocks  (of  magnitude  typically 
m^  = 5)  are  not  simple  events,  but  have  some  tensional 
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subfaulting.  Figure  35  (from  Molnar  et  al.,  1973) 
shows  the  ratios  of  corner  frequencies  from  later  San 
Fernando  earthquakes  ranging  in  magnitude  from  0.5  to 
4.5  Ml. 

Work  with  teleseismic  data  by  Wyss  (1970),  Wyss 
and  Hanks  (1972),  Wyss,  Hanks,  and  Liebermann  (1971), 
Hanks  and  Wyss  (1972)  Wyss  and  Molnar  (1972),  and 
Molnar  and  Wyss  (1972)  has  reached  similar  conclusions 
but  the  Q corrections  are  large  and  uncertain. 

Finally,  Figure  36  displays  the  third  moment  of 
frequency  versus  complexity  discriminant  (from  Anglin, 
1971)  . Anglin  has  indicated  a trend  for  earthquakes 
in  which  increasing  complexity  leads  to  higher  fre- 
quencies. Perhaps,  this  could  be  understood  as  follows 
Complex  earthquakes  imply  many  subearthquakes,  which 
in  turn  implies  a higher  frequency  signal.  A correla- 
tion might  then  also  be  sought  between  complex  earth- 
quakes and  those  close  to  the  limiting  M^m^  line. 

An  understanding  of  short-period  discrimination  is 
crucial  because  despite  several  other  observational 
studies,  e.g.  Booker  and  Mitronovas  (1964),  Lacoss 
(1969),  Ringdal  and  Whitelaw  (1973)  and  Manchee  (1972) 
analysis  mostly  proceeds  in  the  dark  as  far  as  theory 
is  concerned. 
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SUMMARY  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 


We  have  developed  a theory  which  accounts  for  the 
observed  separation  between  Ms:m^  lines  for  earthquakes 
and  explosions,  and  also  explains  the  robustness  of  the 
discriminant  despite  wide  ranges  in  eai  "hquake  (Mg-m^) 
values.  This  wide  range  is  explained  as  being  due  to 
the  existance  of  both  simple  and  complex  earthquakes. 

One  also  sees  that  the  difference  (Mg-m^)  can  be  a 
good  discriminant  even  though  ‘he  slope  of  the  earth- 
quake Mgtm^  line  is  not  1.0. 

The  fact  that  P/S  ratios  are  closer  to  0.3  than 
to  0.03  as  classically  predicted,  together  with  the 
observations  of  higher  P than  S corner  frequency,  have 
been  shown  to  be  related  phenomena  and  explainable  in 
terms  of  tensional  subearthquakes. 

A tentative  explanation  has  been  given  for  the 
success  of  the  complexity- third  moment  of  frequency 
discriminant  (we  might  note  that  this  discriminant  is 
vulnerable  to  evasion  te  ' 'iques) . 

The  results  suggest  .nat  the  projected  intersec- 
tion of  the  earthquake  and  explosion  populations  sug- 
gested by  the  data  of  several  workers,  e.g.  Marshall 
and  Basham  (1972)  , would  not  occur  if  data  at  lower 
magnitudes  were  obtained.  Aside  from  displacement  over- 
shoot, the  only  method  for  obtaining  convergence  of  the 
populations  would  be  for  small  events  to  preferentially 
draw  relaxation  energy  from  small  source  regions  as  in 
the  theory  of  Archambeau  (1968).  This  possibility 
lies  outside  the  scope  of  the  present  tbecry.  Mean 
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^s:mb  re§ressi°n  lines  appear  to  be  somewhat  lacking 
in  physical  significance.  The  fact  that  some  theoreti- 
cal justification  can  be  obtained  for  short-period 
discriminants  suggests  that  they  may  be  relied  upon  in 
non-evasion  situations  if  carefully  regionalized  and  if 
care  is  taken  to  select  paths  unaffected  by  P wave 
multipathing. 

Further  research  should  aim  at  deciding  if  ten- 
sional  subearthquakes  really  exist  or  are  only  a 
theoretical  abstraction.  Field,  experimental,  and 
numerical  work  should  all  have  a place. 

Investigation  of  the  displacement  time  history 
near  the  front  of  the  crack  is  important  to  see  if  on 
the  time  scale  of  interest  telese ismically  (or  regionally 
for  earthquake  damqge  studies)  the  initial  velocity  is 
discontinuous  as  Brune  (1970)  suggested  and  as  we  have 
assumed,,  Brune  (1973)  has  presented  some  experimental 
evidence  to  the  contrary,  but  it  is  not  completely 
clear  that  wave  energy  from  regions  of  earlier  movement 
has  not  slightly  distorted  the  measurements.  This  ques- 
tion is,  or  course,  crucial  because  the  discontinuity  in 
the  displacement  time  history  governs  the  far-field  high- 
frequency  asymptote,  and  any  slope  substantially  different 
from  w would  require  major  changes  in  the  theory. 

The  study  of  the  distribution  of  subearthquakes  is 
also  crucial.  We  have  assumed  a delta  function  here; 
perhaps  N(Mq)  « exp(-MQ)  could  be  shown  to  be  reasonable 
and  tested.  Again,  both  mapping  of  faults  in  the  field 
and  experimental  measurements  are  needed. 
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Although  we  do  not  expect  our  basic  results  to 
change,  it  is  important  to  replace  Haskell's  original 
solution  for  the  simple  earthquake  by  an  exact  solution 
for  a fully  two-dimensional  fault  plane  surface  for  both 
slip  and  tensional  earthquakes.  Here  Savage's  (1974) 
lead  could  be  followed. 


As  a final  point,  Aki  (1967)  and  others  have 
pointed  out  that  for  large  events  m^  is  not  a satisfac- 
tory measure  of  the  total  1„0  Hz  energy  in  the  signal. 

A suitable  measure  suggested  by  Aki  is  the  maximum 
amplitude  in  the  P coda  times  the  "length"  of  the  coda 
raised  to  some  power  between  0.5  and  1.0. 
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Figure  1.  B functions  for  a slip  earthquake  with  slip  in 
direction  of  fault  propagation.  Curves  are  the  same  as 
those  plotted  by  Haskell  (1964) . 
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Figure  5.  Variations  in  amplitude  as  a function  of  propagation 
velocity  for  the  B functions  for  a tensional  earthquake.  Inter- 
mediate slope  is  not  so  apparent  as  in  Figure  3. 
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Figure  8.  S and  P displacement  amplitude  spectra  and  their 
' ratio  for  three  simple  earthquakes  of  different  sizes.  Note 

the  convergence  cf  amplitude  at  high  frequencies;  and  the 
"periodic"  behavior  of  the  ratios. 
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Figure  10.  S and  P displacement  amplitude  spectra  and  their 
ratio  for  three  types  of  simple  earthquake,  plus  one  sublevel. 
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Figure  11.  S and  P displacement  amplitude  spectra  and  their 
ratio  for  a simple  earthquake  plus  0,  1,  2,  5,  or  4 sublevels. 
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Figure  12.  S and  P displacement  amplitude  spectra  and 
their  ratio  for  a simple  earthquake  plus  one  sublevel 
with  a range  of  submoments. 
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Figure  13.  S and  P displacement  amplitude  spectra  and 
their  ratio  for  a simple  earthquake  plus  one  sublevel 
with  a range  of  earthquake  sizes  but  constant  total 
moment . 
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theii^ratio  fn!"!  P.  displacement  amplitude  spectra  and 
tneir  ratio  for  a simple  earthquake,  plus  one  sublevel 
for  a range  propagation  velocities.  J 
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Figure  16.  S and  P displacement  amplitude  spectra  and 
their  ratio  for  a simple  earthquake  and  with  a range  of 
fractional  stress  drops.  An  intermediate  slope  of  u> 
can  be  seen. 
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Figure  18.  20-sccond  total  displacement  versus  1-second 
compress ional  displacement  for  simple  earthquakes  for  a 
range  of  lengths. 


Figure  19.  20-second  total  displacement  versus  1-second 
compress ional  displacement  for  simple  earthquakes  for  a 
range  of  lengths  plus  0,  2,  4 or  6 sublevels  with  constant 
total  moment  leading  to  an  increase  of  25"o.  Ratio  of 
subearthquake  lengths  r = 0.55. 
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•igure  20.  20-second  total  displacement  versus  1-second 
;omp res s ional  displacement  for  simple  earthquakes  for  a 
range  of  lengths,  plus  0,  2,  4 or  6 sublevels  with  constant 
total  moment  leading  to  an  increase  of  1-6.  Ratio  of 
subearthquake  lengths  r = 0.53. 
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Figure  21.  20-second  total  displacement  versus  1-second 
compressional  displacement  for  simple  earthquakes  for  a 
range  of  lengths  plus  0,  2,  4 or  6 sublevels  with  constant 
total  moment  leading  to  an  increase  of  25%.  Ratio  of 
subearthquake  lengths  r = 0.57. 
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Figure  22.  20-second  total  displacement  versus  1-second 
comp  re s s i o n a 1 displacement  for  simple  earthquakes  for  a 
range  of  lengths,  plus  0,  2,  4 or  6 sublevels  with  constant 
total  moment  leading  to  an  increase  of  1%.  Ratio  of 
subearthquake  lengths  r = 0.S7. 
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Figure  25.  20-second  total  displacement  versus  1-second 
congressional  displacement  for  simple  earthquakes  for  a 
range  of  lengths,  plus  4 sublevels  with  a constant  total 
moment  leading  to  an  increase  of  II.  The  effects  of  a 
range  in  fractional  stress  drop  and  propagation  velocity 

can  be  seen. 
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Figure  26.  20-second  total  displacement  versus  1-second 
compress ional  displacement  for  simple  earthquakes  for  a 
range  of  lengths,  plus  4 sublevels  with  constant  total 
moment  leading  to  an  increase  of  II.  The  effects  of 
variation  in  the  type  of  the  simple  earthquake  can  be 
seen  to  be  relatively  minor  except  earthquakes  already 
near  the  limiting  line. 
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Figure  27.  Displacement  amplitude  spectra  for  explosions 
in  granite  using  the  data  of  Werth  and  Herbst  (2963)  and 
the  formulas  of  von  Seggern  and  Blandford  (1972). 
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Figure  23.  20-second  displacement  versus  1-second  displace- 
ment ior  explosions  in  granite  using  the  data  of  Werth  and 
Herb st  (1963)  and  the  formulas  of  von  Seggern  and  Blandford 
(1973),  for  a range  of  yields. 
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Figure  29.  20-second  displacement  versus  1-second  displace 
ment  for  explosions  in  granite  using  the  data  of  Werth  and 
Herbst  (1963)  and  the  formulas  of  von  Seggern  and  Blcndford 
(1973),  for  a range  of  yields.  Figure  has  been  calibrated 
to  Ms  and  values  by  use  of  solution  for  100  kt  explosion 
given  by  Douglas  et  al.  (1973).  Limiting  earthquake  line 
has  also  been  drawn  in. 
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Figure  30.  Ms:mb  pl°t  from  Douglas  et  al.  (1973) (see  Table  2), 
together  with  the  granite  explosion  line  from  Figure  29,  and 
the  limiting  earthquake  line. 


Figure  31.  Ms : m^  for  U.  S.  earthquakes  and  explosions. 
Open  symbols  arc  eartliquakes,  closed  symbols  are  explo- 
sions. From  Evernden  et  al.  (1971).  Granite  explosion 
line  and  limiting  earthquake  line  have  been  drawn  onto 
the  original  figure. 
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Figure  32a.  Pacoima  Dam  accelerogram  of  event  1 
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Figure  32b.  Pacoima  Dam  accelerogram  of  events  2 and  3. 


Figure  32.  Pacoima  Dam  accelerograms  of  San  Fernando 
aftershocks . 
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Figure  32c.  Pacoima  Dam  accelerogram  of  events  4 to  7 
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Figure  3 2 d . Pacoima  Dam  accelerogram  of  events  8 to  10 
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Figure  32f.  Pacoima  Dam  accelerogram  of  events  16  and  17. 
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Figure  32g. 
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Pacoima  Dam  accelerogram  of  events  22  to  25, 
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Figure  ?2h.  Pacoima  Dam  accelerogram  of  events  30  and  31. 
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Figure  34.  P wave  corner  frequencies  versus  S wave  corner 
frequencies  for  large  aftershocks  of  the  1971  San  Fernando 
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earthquake  by  the  strong-motion  instrument  at  Pacoima  Dam. 
From  Molnar  et  al.  (1973)  using  data  from  Trifunac  (1972). 
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Fieure  3b.  Average  third  moment  of  frequency  against 
complexity  for  a population  of  twenty-seven  Eurasian 
explosions  (solid  dots)  and  fifty-one  shallow  Eurasian 
earthquakes  (open  circles).  From  Anglin  (19  1). 


